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INTRODUCTION 

Considerable  effort  has  been  devoted  to  the  solution  of  the 
water-entry  prjb?.em  as  showii  in  reference  1,  v/hich  lists  135 
reports  on  this  subject.  The  d;heories  presented  in  prior  work 
are  in  serious  conflict j  even  for  simple  shapes  such  as  cones. 

For  example.  Hoover  and  Reardon,  reference  2,  use  ':he  assumption 
attributed  to  von  Karman  that  the  amount  of  added  mass  gained 
during  impact  depends  only  upon  the  size  of  the  eody  and  not  at 
all  upon  the  nose  shape,  vihile  Shiffman  and  Spencer,  reference  3^ 
shovj  that  for  c^ne -shaped  bodies,  the  amount  of  added  mass  is  a 
widely  varying  functJon  of  cone  angle.  This  confusion  of  theories 
has  been  compounded  by  c.  relative  lack  of  experimental  data,  v/ith 
the  exception  of  the  work  cf  Watanabe,  reference  4;  Nfejer,  reference 
5j  and  Mosteller,  reference  6.  These  early  (prior  to  I96C) 
experimenters,  however,  were  severely  restricted  by  available 
instrumentation  which  limited  impact  velocity  and/or  configuration, 
and  inadequate  computing  facilities  which  complicated  the  data 
analysis  and  essentially  limited  the  practical  number  of  tests 
conducted.  During  the  last  ten  years,  these  difficulties  have 
to  a  large  extent  been  eliminated.  Furthermore,  facilities  such 
as  the  large  Hydroballistics  Facility  and  associate  Pilot  Tank  at 
the  Naval  Ordnance  Laboratory  have  been  constructed  specifically 
for  the  study  of  water-entry  phenomena.  A  program  was  initiated 
at  NOL  in  1967  to  obtain  basic  vertical  water-impact  acceleration 
data  for  conical  shapes  and  to  reduce  the  results  to  a  simple 
form  for  engineering  applications.  The  conical  shape  was  chosen 
because  of  similitude  considerations.  These  indicated  that  for 
any  one  cone  angle  the  drag  coefficient,  based  upon  the  area 
intersected  by  the  original  water  surfa.ce,  remained  constant;  and 
that  a  linear  relationship  exists  between  the  distance  penetrated 
and  the  effective  distance  penetrated. 

TEST  PROCEDURE 

All  experiments  were  conducted  in  a  similar  manner  with  the 
majority  of  the  tests  being  made  in  the  Pilot  Hydroballistics 
Facility  at  the  Naval  Ordnance  Laboratory  (Fig. 1).  The  model, 
containing  a  single  axially  mounted  crystal  accelerometer,  was 
launched  from  an  air  gun  (velocity,  30  feet  per  second  and  higher) 
or  dropped  through  a  tube,  such  that  normal  impact  with  the  surface 
of  the  water  resulted.  A  cable  (trailing  wire)  connected  the 
accelerometer  within  the  model  to  fixed  electronics.  The  final 
output  of  the  gage  was  photographically  recorded  from  a  CR 
oscilloscope.  As  the  model  neared  the  water  surface,  a  trigger 
screen  was  interrupted,  causing  a  microsecond  strobe  lamp  to 
produce  three  flashes  equally  spaced  in  time,  thereby  exposing 
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the  film  in  an  open-plate  camera.  Located  beside  the  camera 
was  a  photo-pickup  that  converted  the  light  from  the  strobe  into 
voltage  that  v;as  electronically  mixed  with  the  accelerometer  output 
in  the  oscilloscope.  The  arrangement  of  equipment  for  these  tests 
is  shown  in  Figure  2.  The  height  of  the  trigger  was  adjusted 
such  that  the  reflection  of  the  model  on  the  water  surface  was 
photographed,  as  well  as  the  model.  This  technique  resulted  in  tv;o 
photographs,  one  showing  the  model,  generally  above  the  water 
surface,  at  three  points  in  time;  the  other  showing  the  voltage 
output  of  an  accelerometer,  and  the  time  of  each  exposure  on  a 
common  time  base,  displayed  on  an  oscilloscope.  Figure  3  shows  a 
sample  of  these  data.  Additional  data  samples  are  given  in 
Appendix  A. 

The  models  used  in  this  test  series  were  generally  simple 
3-inch  base  diameter  cones  of  various  half  angles  with  cylindrical 
afterbodies  in  which  the  accelerometer  was  mounted.  In  some  cases 
a  sealing  nut  in  the  base  of  the  afterbody  was  employed  to  insure 
a  waterproof  condition  for  the  gage.  The  use  of  the  seal  increased 
the  model  weight  by  about  100  grams.  Certain  of  the  experiments 
involved  models  of  different  construction  when  high  deceleration 
effects  were  desired.  A  special  model  with  a  1.5-inch-diameter 
changeable  nose,  followed  by  a  1.5-inch- diameter  cylinder  and  then 
a  pair  of  3-inch-diaraeter  drag  plates  was  used  for  some  experiments 
involving  slender  cones.  The  high  drag  of  this  model  after  entry 
caused  rapid  deceleration  so  that  it  did  not  become  damaged  when 
it  hit  the  tank  bottom.  This  permitted  multiple  tests  on  a  single 
model  and  reduced  the  overall  cost  of  the  program.  A  20-degree 
aluminxim  nose  on  this  model  showed  no  damage  when  launched  at 
velocities  greater  than  500  feet  per  second  into  less  than  8  feet 
of  water.  In  all,  twenty  different  models  were  tested.  The 
model  details  are  shown  in  Appendix  B. 

Three  accelerometers  were  used  during  the  test  series: 
a  high-capacitance,  high-impedance  multiple -crystal  gage;  and 
two  similar  quartz  gages  possessing  internal  electronics  and  low 
output  impedance.  The  high  capacitance  of  the  multiple-crystal 
gage  reduced  the  relative  error  caused  by  changes  in  capacitance 
in  the  cable  (trailing  wire)  due  to  acceleration  and  other  loads. 

Quartz  gages  were  used  in  the  lightweight  models,  as  well  as 
some  of  the  lighter  than  normal  heavy  models.  An  electronic 
low-pass  filter  was  often  used  with  these  gages  to  eliminate 
ringing  of  the  transducer  or  model. 

The  large  capacitor  crystal  accelerometer  was  used  with  a 
charge  amplifier  to  increase  the  system  RC  time  constant  when  a 
maximum  acceleration  of  less  than  30g  was  expected.  An 
electronic  filter  was  also  frequently  used  with  this  gage. 

The  gages,  scopes,  filters,  and  charge  amplifier  were 
calibrated  as  a  system.  Several  combinations  were  recalibrated 
during  the  test  series.  Three  methods  of  calibration  were  used 
for  each  gage  and  electronics  combination  employed  in  the  test 
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series.  The  three  calibration  methods  vrere:  (1)  shaker  table 
with  optical  measurement  of  displacement;  (2)  shaker  table  with 
comparison  with  a  standard  gage;  and  ^3)  a  one  g  drop. 

In  the  one-g  drop  test,  the  ^ge  and  model  were  held  by  a 
short  lengoh  of  nichrome  wire.  A  high  current  electric  source 
was  applied  to  the  wire,  causing  rapid  melting  and  quick  release. 
This  step  pulse  of  one  g  was  used  to  determine  the  leakage  rate 
of  the  gage-electronic  system,  as  well  as  to  determine  the  gage 
constant. 

The  gage  constants  used  for  each  system  were  the  simple 
averages  of  the  values  determined  from  all  of  the  calibration 
techniques . 


RANGE  OF  VARLABI£S  TESTED 


In  the  test  series,  a  wide  range  of  variables  was  investigated 
experimentally  to  establish  theoretical  postulations  such  as 
similitude  and  added  mass  concepts.  In  general,  3-inch-diameter 
models  were  tested,  except  for  some  1.5-inch-diameter  noses  as 
explained  in  the  previous  section.  This  limitation  in  maximum 
diameter  resulted  from  the  maximum  barrel  diameter  available  for 
the  air  gun  used.  It  is  hoped  that  future  tests  v;ill  be  conducted 
on  much  larger  diameter  models  to  investigate  the  conclusion 
reached  herein  that  similitude  does  exist  for  cones.  The  ranges 
of  test  variables  v;ere  as  follows: 

1.  total  cone  angle  -  10  to  140  degrees 

2.  velocity  -  10. 5  feet  per  second  to  250  feet  per  second 

3.  mass  -  128  grams  to  9387  grams. 

The  distribution  of  these  variations  is  given  numerically  in 
Appendix  C  and  graphically  in  Figure  4. 

DATA  REDUCTION 

The  "a  priori"  assumptions  made  in  the  data  reduction  were 
that  the  external  forces  depend  only  upon  size,  geometry,  actual 
velocity,  and  distance  below  the  water  surface;  and  that  the  law 
of  similitude  holds  for  cones.  Additionally,  it  was  assiuned  that: 

1.  buoyancy  force  is  calculated  from  original  water  surface; 

2.  friction  drag  coefficient  is  constant  for  all  conditions; 

3.  steady-state  drag-coefficient  is  constant  for  each  cone  angle 
and  equal  to  the  value  obtained  in  water-tunnel  tests  extrapolated 

to  zero  cavitation  number  (ref.  7  for  0-45  degrees,  ref.  8  for 
60-l40  degrees.  Table  1);  and 
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4.  surface  heave-distance  is  proportional  to  the  depth 
of  penetration  for  each  cone  angle. 


Table  1 

Steady-State  Drag  Coefficient 


Cone 

.4ngle 

Helena rdt 
Ref.  (8) 

Hoerner 
Ref.  (7) 

Value 

Used 

10 

.05^ 

.0595 

15 

.0684 

.0584 

20 

.0866 

.0866 

30 

.152 

.1410 

.141 

45 

.270 

.257 

.257 

60 

-365 

.^15 

-365 

90 

.497 

.^97 

120 

•  596 

.596 

140 

.660 

.660 

180 

.79 

{-T9) 

The  preliminary  data  reduction  was  done  manually  and  involved 
changing  the  photographic  records  to  digital  form.  The  photo¬ 
graphic  information  on  the  plate  camera  and  oscilloscope  records 
was  reduced  to  numbers  as  follows:  the  plate  camera  photograph 
was  used  to  determine  the  impact  velocity  and  to  locate  the  water 
surface  relative  to  the  oscilloscope  trace.  The  distance  the 
model  advanced  between  flashes  was  computed  using  the  diameter 
of  the  model  as  a  reference  dimension.  The  time  between  flashes 
was  measured  by  an  electronic  counter  and  the  distance  from  the 
reference  picture  was  measured  (if  possible)  to  a  point  half-way 
between  the  picture  and  the  matching  reflection.  Here  again, 
the  model  diameter  -was  used  as  the  reference  length. 

The  readings  of  the  oscilloscope  records  were  made  using 
a  film  reader  or  toolmaker's  microscope  in  the  following  order: 
the  distance  of  a  vertical  cm  was  read  to  determine  the  scale 
factor j  next,  the  pip  relating  the  model's  position  v?ith  the 
water  surface  v/as  readj  and  then  a  series  of  points  along  the 
trace  were  read.  Both  components  of  each  point  were  recorded. 

The  numbers  obtained  from  the  pictures,  along  with  recorded 
information  such  as  model  diameter,  weight,  and  scope  gain, 
were  entered  into  a  time-shared  digital  computer. 

The  acceleration  time  data  were  converted  to  English 
engineering  units  for  computations;  and  when  needed,  a  scope 
distoration  was  applied.  The  coordinate  axis  was  shifted  such 


8 


NOLTR  71-25 


that  zero  acceleration  occurred  just  before  impact.  The  follov;ing 
correction  for  leakage  was  then  applied: 


dt 


(True) 


dU 

rF 


(read) 


RC  Jq 


(read)°^ 


The  coordinate  axes  were  then  shifted  such  that  -f-  ig  occurred 
just  above  the  water  and  that  time  zero  occurred  at  water  contact. 
The  velocity  and  distance  were  then  determined  for  each  data 
point  using  trapezoidal  integration.  The  values  related  to  v:ater- 
entiy  effects  were  then  computed. 


The  equation  relating  the  ssomentum  just  before  first  contact 
v;ith  the  momentum  at  some  later  time  is  given  as  Equation  (1), 
and  the  force  equation  is  given  as  Equation  (2); 


n  M-U(!te)  =  B  dt  -Mgt  +  IJJ  AU®dt 


-  H  (««■)  -  U  I?  =  3-Mg  I  C,,  AU® 


2  ds 


(1) 

(2) 


Substituting  ^  for  the  second  term  of  Equation  (2)  and 
rearranging  tertas'^gives : 

dU  r,  ^  _ dm  ,,2  P  -  2 

-  St  -  B  4-  —  V  2  C^g  Av 


(3) 


If  the  total  drag  coefficient  is  defined  as 

„  _  2  dm  .  ^ 

'"d  oA  dS  '  ^ds 


(^) 


then  Equation  (3)  inay  be  rewritten  as 

-  ^  (?4-Hn)  -  B  +  J4g  =  I  D  A'/  (5) 

The  value  of  added  mass  (m)  and  ^  can  be  calculated  using 
Equation  (3)  and  the  total  drag  coefficient  v/as  computed 
from  Equation  (4).  Hov/ever,  a  simpler  procedure  is  to  determine 
m  from  (1)  and  C.  from  (5).  In  order  to  ensure  an  accurate 
evaluation  of  the  integrals  in  Equation  (1),  the  added  mass  (m) 
v;as  also  computed  by  integrating  ^  as  given  in  Equation  (4) 
using  the  computed  value  of  C^. 


For  the  vertical  entry  of  cones,  the  added  mass  should  be 
given  by 


m  = 


c;  -a  o 

ok  (tan 


(6) 


Then 
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g  =  3p5c  (tan  (7) 

Substituting  (7)  into  (3)  yields 

-  if  -  B  +  Mg  =  (3olc(tan  |)  V  -f  |  C^^A)\J^  (8) 

If  the  total  mass  constant  (K)  is  defined  as: 

K  =  k  +  (tan  ^  (9) 

then  (8)  can  be  rewritten  as 

-  ^  (M-Ha)  -  B  +  ^ig  =  3PK(tan  S^U^  (10) 

This  form  of  the  equation  was  chosen  to  compare  easily  the 
experimental  results  with  the  theoretical  results  obtained  in 
reference  3* 

It  should  be  noted  that  the  total  mass  constant  (K)  is  also 
a  drag  coefficient  related  to  the  area  intersected  by  the  original 
water  surface.  The  relationship  between  the  mass  constant  (K)  and 
the  total  drag  coefficient  Ca  is  obtained  by  subtracting  Equation 
(5)  from  Equation  (10)  and  dividing  out  p  and  U^: 

=  6K(tan  |)3s^(A)"^  (11) 

If  the  area  used  to  relate  to  drag  force  is  the  area  intersected 
by  the  heaved  water  surface, “then  the  instantaneous  drag  coefficient 
is  given  by 

=  (6K  tan  |-}{tt)  ^  (12) 

It  was  expected  that  K,  and  hence  would  be  a  constant  during 

entry,  as  suggested  in  reference  3*  investigate  this,  a  straight 
line  least  square  approximation  was  fitted  to  the  computed  values 
of  K  and  S.  The  range  in  S  extended  from  20  percent  to  90  percent 
of  the  depth,  at  which  the  maximum  total  drag  coefficient  (C^j) 
occurred.  The  exclusion  of  the  first  20  percent  eliminated  most 
of  the  points  with  high  relative  reading  errors  and  the  last  10 
percent  was  not  used  because  of  an  observed  "rolling  off"  as  the 
point  of  separation  approached. 

It  is  seen  from  Equation  (10)  that  the  slope  of  the  mass 
constant  (K)  versus  depth  (S)  curve  can  be  altered  by  changing 
the  value  of  the  depth  (S)  associated  with  each  data  point. 
Variations  changes  in  S  can  be  accomplished  by  adjusting  the 
value  of  the  distance  to  the  water  obtained  from  the  plate  camera 
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FIG .  6  MEASURED  DEPTH  COMPARED  TO  COMPUTER  DEPTH 
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I  film.  This  value  was  systematically  changed  by  the  computer  until 

'  a  movement  of  .0005  o£  the  cone  length  \vould  cause  a  sign  reversal 

I  in  slope.  This  slope  was  then  recorded,  as  was  the  sura  of  the 

f  squares  of  the  residuals  between  the  least  square  fit  and  the 

[I  data  points. 


k 

t 

L' 
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The  computed  depth  was  that  at  which  the  maximum  total  drag 
coefficient  occurred  after  the  least  square  perturbations.  A 
measured  depth  was  obtained  for  many  tests  either  by  comparing 
the  distance  to  the  water  measured  from  the  film,  with  the  value 
arrived  at  by  the  computer,  as  shoxim  in  Figure  6;  or,  in  a  fev/ 
cases,  directly  from  the  photograph  and  scope  trace. 
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The  penetration  ratio  h  was  computed  by  dividing  the 
computed  depth  by  the  cone  length.  A  related  terra  is  the 
surface  heave  defined  as  the  cone  length  minus  the  computed 
depth.  The  numerical  results  for  each  test  are  given  in  Appendix 

C. 


RESULTS  AND  CONCLUSIONS 

Perhaps  the  most  important  result  of  the  entire  test  series 
is  proof  that  the  total  added  mass  constant  K  is  indeed  constant 
as  the  cone  is  wetted.  A  constant  K  implies  from  Equation  (12) 
a  constant  instantaneous  drag  coefficient  Cdj|^,  and  from  Equation 
(9)  a  constant  added  mass  constant  (m).  After  the  perturbation 
in  distance,  the  average  measured  point  was  within  3*6  percent 
(root  mean  square)  of  a  constant  value.  Also,  a  comparison  of 
all  measured  depth  and  corresponding  computer  depths  shovjs  an 
average  difference  of  .5^  percent.  To  the  satisfaction  of  the 
author,  these  results  show  that  K,  m  and  C(j.  are  constants  over 
the  range  of  depths  used  in  the  least  square  fit.  There  seems  to 
be  a  skew  of  about  .3  percent  betv/een  computed  depth  and  measured 
depth,  as  shown  in  Table  2. 


Table  2 

Comparison  of  Measured  Depth  and  Computer  Depth 


Rank 

Cone 

Angle 

Shot  No. 

Difference 

(inches/inch) 

Aggregate 
Error  % 

0 

- 

- 

- 

.54 

1 

30 

1895 

.1467 

.28 

2 

45 

1662 

•  1295 

08 

3 

20 

1902 

-.0927 

.23 

4 

90 

1695 

-.0844 

.40 

5 

120 

1997 

.0822 

.25 

30 

1888 

-.0754 

.40 

7 

i4c 

2024 

-.0735 

.56 
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Table  2  (Cont.) 


Rank 

Cone 

Angle 

Shot  No. 

Difference 
(inch ) 

Aggregate 

Error 

8 

140 

201405 

.0714 

.43 

9-14 

- 

6  pts 

• 

0 

CTi 

1 

» 

0 

-.11 

15-22 

- 

8  pts 

.05-. 06 

+  .49 

The  table  shows  the  effect  of  discarding  data  points  in  descending 
order  of  percent  difference.  Perhaps  this  difference  is  due  to 
a  systematic  reading  error. 

It  is  hoped  that  one  of  the  major  contributions  of  this 
report  is  the  determination  of  total  drag  coefficient  (C^)  for 
most  cone  angles.  To  be  of  maximum  utility,  should  be 
independent  with  respect  to  velocity,  size,  ana  vjeight.  In 
order  to  investigate  the  velocity  dependence,  all  reported  values 
of  (Appendix  C)  were  normalized  and  a  least  square  fit  to 
C(i  =  A+BUq  performed.  This  calculation  resulted  in  =  I.008  - 

.00012Uq.  The  results  should  be  useful  in  sizes  below  and  above 
those  tested  because  of  similitude.  The  author  could  see  no 
systematic  effect  on  due  to  varying  weights.  As  seen  in 
Table  3^  the  drag  coefficient  for  the  90-degree,  3-inch-diameter 
cones  remained  within  2  percent  of  the  average  for  a  model  weight 
range  from  .33  pound  to  20.9  pounds. 

Table  3 

Maximum  Drag  Coefficients  for  Various  Model  Weights 


Cone 

Normal 

— 

Odd  Wt  Cd 

Angle 

VJeight 

Weight 

Av  Wt  C(i 

140 

1900  (6) 

251 

(3) 

.971 

120 

1700  (26) 

128 

.954 

287 

b) 

.949 

90 

1600  (17) 

150 

(3) 

.992 

284 

(2) 

1.011 

9489 

(^) 

1.018 

60 

1800  (8) 

(5) 

.985 

45 

1900  (10) 

(4) 

1.057 

20 

2200  (9) 

295 

(1) 

.935 

No  odd 

weights  were 

tested 

at 

cone  angles 

of  30, 

15.  or  10  degrees 

• 

.  .... 
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Tr.e  c  lur.*  c^r.es  'i2~  ie-^rerrs  ar.d  l^^-C  degrees)  did  indeed  show 
a  rei-^-ior.  of  atc”0  5  peroer.t  for  the  light  aadeis.  This  is 
lelievei  to  ce  due  tc  er’^crs  ir.  the  ^ge  constant  (large  percent 
s-owdrifr.',  zo-del  weight  'how  orueh  wire  to  include),  and  trace 
reading  errors-  The  average  oaxiooio:  total  drag  coefficient  for 
eaor.  oor.e  angle  as  a  funotion  of  this  angle  is  shown  in  Figure  7. 

Also  snown  in  Figure  are  *.*alues  deri\'ed  fros  reference  5  snd 
tnose  ass'uned  for  T-c-,  Table  1.  liote  tne  basic  agreement  of 
va.ues  'f  zaxinun  drag  ooeff’cient  and  the  convergence  of  to 
:'zr  “re  snaller  cone  angles.  The  curve  shows  that  the  totai  drag 
'refflolent  is  approxinateljr  equal  to  the  steady-state  drag 
oceffioient  f roco  aero  cone  angle  to  a  total  cone  angle  of  about 
degrees,  above  which  rapid  divergence  occurs.  Another  v;ay  of 
stating  this  is  that  the  naxirous  shock  factor  during  impact  as  a 
funooion  of  cone  angle  is  almost  unity  up  to  45  degrees,  and  then 
.'apioly  Inoreases  to  values  of  2  and  5  blunter  cones.  Also, 
note  the  excellent  agreement  with  the  theoretical  values  (Wagner 
arjd  fonneider,  reference  2nd  the  general  agreement  with  the 
e'O er iments  of  Majer  and  iiatanabe  (reference  4). 

An  important  result  of  these  tests  is  the  percent  of  the 
cone  surcierged  'cased  upon  original  water  surface)  when  the 
oaxioiom  total  drag  coefficient  occurs.  These  data  are  shovm 
in  Figure  3.  A^in,  reasortable  agreement  occurs  when  compared 
with  previously  selected  tests. 

It  is  seen  that  if  the  cone  angle  is  very  small,  then  the 
2!axis:jci  drag  coefficient  occurs  after  the  cone  base  passes  the 
original  water  surface.  The  number  of  tests  conducted  on  10- 
degree  and  15-degree  total  angle  cones  leaves  this  point  somewhat 
in  doubt.  For  clunter  cones,  the  percent  of  cone  length  penetrated 
drops  as  the  cone  angle  increases  until  a  value  of  71.5  is  reached 
for  1-0 -degree  cones.  This  is  equivalent  to  51*1  percent  of  area 
celow  the  origirtal  water  surface.*  When  compared  tc  the  simple 
theory  —sed  in  reference  2,  this  change  in  radius  couid  change 
the  added  nass  by  almost  three  to  one. 

Oomparison  of  the  test  results  with  theory  is  best  shown 
in  Figure  9,  where  the  experimental  values  of  total  mass  constant 
rl)  are  corrected  for  friction  and  compared  with  the  theoretical 
V® lues  arrived  at  in  reference  3*  Also  plotted  in  Figure  9  sre  ex— 
oerimer.tal  values  cotaxned  from  reference  3.  For  small  total  cone 
angles,  the  cheoretical  values  are  much  greater  than  those  obtained 
from  the  experiments  even  with  the  steady-state  drag  term  added 
to  the  added  mass.  For  very  blunt  cones,  the  theoretical  predic¬ 
tions  seem  to  be  roo  small.  The  values  of  added  mass  constant  derived 
in  reference  3  from  the  tests  described  in  reference  4  are 
monctonic,  decreasing  v;ith  respect  to  increasing  depth.  These 
could  procably  be  brought  into  agreement  v/ith  the  lav/  of  similitude 
by  a  slignt  change  in  time  of  impact. 
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FIG.  7  MAXIMUM  TOTAL  DRAG  COEFFICIENT  VS  TOTAL  CONE  ANGLE 
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The  uncorrected  mass  constants  versus  cone  angle  for  fine 
cones  are  shown  in  Figure  10  as  is  the  effect  of  assuming  a 
friction  drag  coefficient  (C^)  of  .OO3I.  It  is  seen  that  the 
subtraction  of  skin  friction* causes  the  curve  of  the  mass  constant 
to  pass  through  zero  for  zero  cone  angle.  The  friction  effect 
was  also  used  to  correct  the  total  drag  coefficient  C^,  see 
Figure  7. 

In  general,  for  each  cone  angle,  C(j,  k,  K,  and  m  did  not  seem 
to  vary  with  velocity  or  mass  of  the  model  over  the  range  of 
variables  tested. 

Another  important  result  is  the  amount  of  added  mass  accumu¬ 
lated  during  entry.  The  experimental  values  of  added  mass  (m) 
associated  with  maximum  drag  coefficient  versus  cone  angle  are 
shewn  in  Fig>ire  11.  As  shovm  in  Figure  5>  the  final  results 
regarding  added  mass  are  strongly  dependent  on  the  assumptions 
made  concerning  C^jg.  Changes  in  these  asstunptions,  if  desired  in 
the  future,  may  be  made  analytically  as  seen  from  Figure  5* 

Note  that  the  added  mass  constant  (k)  v;ill  not  give  the  same 
exact  values  of  ra  due  to  the  roll-off  experienced  as  the  point 
of  maximum  total  drag  coefficient  is  approached.  The  added  mass 
constants  defined  by  the  equations: 

m  =  DkS^  (13) 

and 

m  =  ok(|-)^  (1^) 

are  shown  versus  total  cone  angle  in  Figure  12,  The  average 
added  mass  associated  with  the  maximum  total  drag  coefficient, 
the  cone  length,  and  the  average  penetration  ratio  were  used  to 
compute  each  plotted  point - 

As  a  body  enters  the  water,  the  amount  of  added  mass  should 
increase  after  the  maximum  drag  coefficient  is  experienced  due  to 
water  surface  effects.  Tiis  is  shown  graphically  in  Figure  5* 

The  values  of  added  mass  for  points  near  and  after  C(j  for  many 
tests  are  given  in  Appendix  C.  The  added  mass  is  showft  ever  the 
entire  range  of  depth  I’or  three  tests  in  Figure  I3,  These  results 
are  in  qualitative  agreement  with  reference  9*  many  of  the 

tests  with  slender  cones,  the  added  mass  curve  did  not  become 
level  at  as  great  a  depth  as  those  shown  on  Figure  13i  Appendix  C. 
It  is  believed  that  this  prem.ature  reduction  in  added  mass  was 
caused  by  using  unstable  models  that  developed  significant  angles 
of  yaw  after  entry  and  thus  misaligned  the  accelerometer  axis 
with  respect  to  the  velocity  vector. 

The  photographs  taken  by  the  plate  camera  sometimes  contained 
an  image  of  the  model  entering  the  water  in  addition  to  the 
required  above-water  information.  A  study  of  these  splash 
pictures  and  associated  acceleration  traces  leads  to  the  conclusion 
that  the  maximum  deceleration  occurs  as  soon  as  the  base  area 


17 


71-25 


FIG.  10  CORRECTION  OF  TOTAL  ADDED  MASS  CONSTANT 
VS  TOTAL  CONE  ANGLE 


FiG .  1 1  ADDED  MASS  ASSOCIATED  WITH  MAXIMUM  DRAG 
COEFFICIENT  VS  CONE  ANGLE 
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FIG.  12  ADDED  MASS  CONSTANT  VS  TOTAL  CONE  ANGLE 
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bears  u^on  solid  water  for  the  bltanter  cones  and  somev;hat  later 
for  the  slender  cones.  The  data  from  shots  18^7  and  1687  included 
in  Appendix  A  v;ere  c  'ong  those  used  to  reach  this  conclusion. 

APPLICATION  OF  RESULTS 

The  practical  application  of  these  results  will  yield  the 
pvedicclon  -^f  the  axial  forces  experienced  by  conical  bodies 
during  vertical  water  entry.  In  computing  these  forces.  Equation 
(15)  nay  be  ased: 

-  II  (2.1-rm)  -  B  +  Mg  =  I  DCd(S)A.pU^  (15) 


or  in  another  form: 


(16) 


Note  that  the  internal  body  forces  and  equal  external  pressure 
forces  (axial)  are  given  by: 


2^  +  JP  dA  =  P  =  M  II 


(17) 


For  any  particular  design  known  constants  include:  total 
cone  angle  (0),  missile  mass  (M|,  base  radius  (R),  and  initial 
velocity  (Ug).  In  Equation  (I6),  the  remaining  functions  are: 
added  mass  (m),  buoyancy  (B)  and  total  drag  coefficient  (C^(S)). 

Due  to  the  shape  of  these  functions,  it  is  necessary  to  determine 
their  values  over  two  ranges  of  depths.  From  contact  until  the 
depth  at  v/liich  the  total  drag  coefficient  is  obtained.  The  functions 
for  added  mass  m  and  total  drag  coefficient  are  given  by 


tn  =  pkS' 

C<i(S)  =  {Cd 


max 


4)") 


For  depths  greater  than  Lh,  these  functions  are  given  by 

=  Plc{Lh)3  +  ^ 


m 


A 


Cd^S)  =((Cd  )  -  Cds)  173 

/VJO  -v  — ' 


max 


'ds 


where 


A  =  exp[-  2L3§imzM!l-i  +  .3.ri  .  (^) 

Cj(S)  in  this  case  is  a  curve  fit  to  the  data.  The  values  of 
^ds  ^sed  in  Equations  (19) ^  (20),  and  (21)  should  i 


(18) 

(19) 

(20) 
(21) 

(22) 

nclude 
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the  effect  of  skin  friction.  For  example,  if  the  Reynolds 
number  is  about  6  x  10°,  then  .OO3  cotan  ^  should  be  added  to  the 
values  obtained  from  Figure  7- 

The  values  of  k  used  in  Equations  (I8)  and  (20)  are  obtained 
from  the  upper  curve  of  Figure  12. 

From  contact  until  the  cone  base  is  level  with  the  original 
water  surface  the  buoyancy  is  given  by 

B  =  (tan  Ifs'i  (S3) 

and  for  greater  depths  (S>L) 

B  =  .DgrR2(s  -  I  L)  (24) 

Note  these  equations  neglect  "roll  off",  which  should  give 
negligible  errors.  If  greater  accuracy  is  desired,  use  m 
derived  from  K  up  to  80  percent  immersion,  and  a  smooth  function 
to  the  peak  value.  The  only  remaining  unknovms  in  Equation  (16) 
are  U  and  ^  which  can  be  solved  for  by  various  numerical  methods . 
One  of  the^Simplest  assumes  that  the  velocity  U  is  for  the 
first  try,  and  U  (from  a  preceding  try)  for  each  succeeding 
iteration.  When  U  is  determined,  then  the  acceleration  ^  may  be 
solved  for  directly  using  Equation  (5).  A  simple  computer  program 
for  solving  this  problem  is  given  in  Appendix  D.  Example  runs 
for  a  3-ihch,  5-potind,  90-degree  cone  and  an  l8-inch,  l8-pound, 
90-degree  cone  are  also  included.  The  relatively  large  decelera¬ 
tions  of  the  large  light  cone  caused  the  maximum  deceleration  to 
occur  before  the  maximum  total  drag  coefficient. 

RECOMMENDATIONS 

Based  upon  the  experience  gained  from  this  test  series  and 
the  results  thereof,  the  author  recommends  in  addition  to 
confirming  tests,  that  the  following  work  be  conducted  regarding 
vertical  water  entry  of  cones: 

1.  Axial  forces  experienced  by  cones  during  water  entry 
be  extended  to  significantly  higher  velocities . 

2.  Assumptions  of  buoyancy  and  C^s  confirmed:  buoyancy, 
by  very  large  models;  and  C^jg,  by  a  method  suggested  in  reference 
10,  by  testing  two  geometrically  similar  models  of  differing  mass, 
and  solving  simultaneously  for  and  m  using  Equation  (3). 


3.  Effects  of  angle  of  attack  at  entry  be  determined. 

4.  Investigation  as  to  why  the  added  mass  maximizes  before 
complete  cavity  development  be  made. 
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APPENDIX  A 


? 

I 


d 


The  original  data  from  several  tests  form  this  appendix. 
These  particular  tests  were  selected  to  provide  information  at 
most  cone  angles  for  both  heavy  and  light  models.  Shot  numbers, 
weight,  gage  constants  and  other  necessary  numeric  values  are 
given  in  Table  A-1.  The  remainder  of  this  appendix  contains 
optical  copies  of  the  original  data  photographs. 
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APPENDIX  B 


This  section  consists  of  sketches  depicting  the  various 
models  tested  during  this  series.  The  measured  variation  from 
the  nominal  values  used  to  reduce  rhe  data  vjere: 

3-inch  base  diameter  +  .000  -  .OO3  inch 
1.5-inch  base  diameter  +  .001  inch 

Cone  angle  +  .1  degree  for  cone  angles  of  10  through  30  degrees 
and  l40  degrees 

+  .3  degree  for  cone  angles  of  ^3  through  120  degrees 

As  the  test  series  progressed,  a  slight  rounding  of  some  cone  tips 
occurred,  due  to  repeated  hitting  of  the  nylon  impact  mats  used  to 
stop  the  models.  The  maximum  tip  radius  was  measured  as  .O35  inch. 
As_ mentioned  before,  the  use  of  the  drag  plate  model  (Figure  B-2) 
eliminated  noticeable  tip  damage  for  aluminiim  cones  as  fine  as 
10-degree  total  angle.  No  correction  was  attempted  for  the 
rounding  of  the  cone  tips. 
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APPENDIX  C 


In  this  appendix  are  numerical  values  obtained  for  each  test 
in  this  series.  Also  included  are  plots  of  average  values  of 
total  drag  coefficient  and  added  mass  versus  depths  for  several 
cone  angles. 


The  following  terms  have  not  been  previously  defined,  at  least 
as  used  in  these  tables. 

o 

added  mass  coefficient  k;  in  m  =  pk(S/h) 
surface  heave  -  cone  length  minus  computer  depth 


residuals  -  sum  of  the  square  of  the  residuals  resulting  from 
the  least  sqiiare  fit  of  mass  constant  (total).  The 
numbfer  in  parentheses  is  the  number  of  data  points  used 
in  the  least  square  fit.  The  average  residual  is  the 
square  of  the  arithmetic  average.  The  percent  residual 
is  100  (average  residual/average  total  mass  constant). 


Cd(f)  -  the  amount  of  Cd  delegated  to  skin  friction 
ca  ratio  - 

Cd 


In  those  cases  where  the  tabulated  number  is  not  dimensionless, 
the  property  listed  is  for  a  3-inch-diameter  body  except  for  10- 
and  15-degree  cones  where  the  values  listed  appoint  a  shot  number 
relating  to  1.5-inch  diameter. 
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Table  C-l.  SUMMARY  OF  CONE  DATA 
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Table  C-3.  20-DEGREE  TOTAL-ANGLE  CONES 
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FIG.  C-5  AVERAGE  TOTAL  DRAG  COEFFICIENT  AND  ADDED  MASS  vs  DEPTH  FOR  140°  CONES 
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APPEiroiX  D 


This  computer  program  written  in  BASIC  provides  the  nuicerioal 
solution  to  the  force  equation  associated  v;ith  the  v?ater  entiy  of 
cones.  For  the  sake  of  simplicity,  the  drag  coefficients  and  the 
penetration  ratios  are  required  inputs  for  each  run,  as  is  the 
cone  angle.  If  desired,  Cdjj^axi  Cd*s  and  h  may  be  easily  computed 
internally  from  polynomial  equations  in  8.  The  various  sections 
of  the  program  perform  the  function  as  follov;s: 

SUB  500  -  reads  in  data  and  establishes  the  depth  stations 

SUB  700  -  computes  buoyancy  and  total  drag  coefficient  for 
each  depth  station 

SUB  900  -  computes  added  mass  for  each  depth  station 

1000  -  computes  dU/dS  based  upon  previous  velocity  and 
integrates  to  obtain  new  velocity.  The  interaction 
process  is  stopped  v/hen  is  sufficiently  small 

SUB  1200  -  computes  acceleration 

SUB  1500  -  prints  out 

Line  I700  -  data  input 

The  data  input,  in  order,  consists  of:  maximum  total  drag 
coefficient  (Cd);  steady-state  drag  coefficient  (Cds)^  added 
mass  constant  (k)^  and  total  cone  angle  (9)  in  degrees;  base 
diameter  in  feet;  penetration  ratio  (h);  body  weight  in  pounds; 
and  entry  velocity  (Uo)  in  feet  per  second. 

NOTE:  Do  not  forget  to  include  skin  friction  in  maximum  total 
drag  coefficient  and  steady-state  drag  coefficient. 
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no  GO  S'JS  700 

IIS  GO  S'Jo  900 

120  GO  S«ja  1000 

12s  G3  S’Jd  1200 

130  G3  SUB  1500 

500  head  C1,C2^K1,AA 

505  HEAD  D2>H,Wl 

507  HEhD  91 

510  LET  V2=0 

515  LET  V3=0 

520  LET  SCI )=0 

5-ii5  LET  A4=Aa/C2^'57.3) 

550  LET  A5=.5^'D2/TANCA4) 

554  LET  A5=A5^=r{ 

555  LET  D1=A6/15 
555  E9H  r<=2  TO  30 

560  LET  SCK)=SCK- 1 )+Dl 
56?  !\‘E\T  K 

565  LEI  SC31 )=SC30)+5*D1 
570  LET  SC32)=SC30)+10»D: 

575  LET  SC33)=S(30)+2C=!'01 

5-10  HETUhN 

700  LET  SC !)=0 

705  LET  CC 1 )=0 

710  FOR  K=1  TO  33 

715  LET  HCK)=62.4-:C3.1416/3)*SCK)t3=»=CTANCA4))t2 
720  IF  SCK)<AS  GO  TO  730 

725LET  3CK)=62. 4»C3.  I  416/4>*D2t2»CSCK)-C2/3)=fA5) 

730  NEXT  K 

735  FOR  K=1  TO  16 

740  LET  CCK)=C1*SCK) T2/A6t2 

745  NEXT  K 

7  50  LET  H=.  5=^02 

755  F0H  K=17  TO  ?3 

76?  LET  C4=.661*CCSC.K)-SC  16))/SC  16)) 

765  LET  A9=.3/C1+CSCK)-SC  16))/C.5?=D2))»?+EXPC-9*C4) 

770  LET  CCK)=CC1-C2)*C 1/1 .3)*A9+C2 

780  NEXT  K 

785  FOR  K  =  1  T0  33 

7^^7  LET  VCK)=V1 

789  NEXT  K 

790  kETUHN 
900  LET  MC 1 )=0 
910  FOR  K=1  TO  16 

920  LET  MCK)=C62.4/32.2)=fK1tSCK)  t3 

9  30  NEXT  K 

940  ~0R  K=17  TO  33 

945  LET  P=SCK)-SCK-1) 
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V5U  Ldl  &  In  )=t‘lK  -  I  J-H  Cir.  A/  (  A.  Ay>  .  lb  bA*  bi-T  ir'  C  .  <  C  IK  >  *Ct  K  -  1  ^  ^ 

9  f.u  eVLaI  .k 
9  7U  hclUHN 
1000  FOn  K=?  10  33 

1010  Lcl  ZIK.)=-  (.(.k.  A/{.A.  AJ^HKi^  .IbbA^Lbib^’  Kii^.y 
101b  Lcl  =  <C<K)  +  Cr  <K1-Un/ VtK  >  >/(  t  LI /3ir.ir  z+MK.- . 

NcAl  K 

1030  FOn  K=  i-lO  33 

lOAO  Lcl  V<K)  =  V  IK- 1  5^  <4  (K) +ZtK- 1  >  >>  C  i(K  >-LCrt- 1  y  ^ 

lubO  NcAl  K 

10^.0  Lcl  Vb=Mi  1  ( Vi-- V<  1  r  >  i+rtLlt  V3- VC  3^- ) ) 

1070  IF  V3<.001  IhLN'  lllu 
lObO  LLl  Vi:=VC  1  Ci 
1090  LLl  V  3=  VC  Si”) 

1:100  fc  1C  lUOO 

1110  hblUhW 

IfeOO  FCK  K  =1  10  33 

iiflO  Lcl  k(.H:  =  -  I  C  (k.  A/  {  A.  A^=^  .  7b  b4=^C  C  K  )  »  Lbt  V- CK  >  t  b  +  L  C  K  i  -  H  / 

Iblb  LLl  «CK)=hCK)/cc^  l/3b.b)+J' <K>) 

IbbO  NcAl  K 

1  b30  KclUKN 

1500  Lcl  hb  =  k*bl.Z>t^A 

1505  f"Klr-.  l  '  COiVL  ANPLi:.  =  *  nb 

1511/  HKIMI  •LilA.t^clcn  IW  FI  =’ tb 

I5b0  PKIWI  Lcirhl  hV  Ltl  =•  I  1 

1530  PKIXI  ■  L-lil*  >•  VcL’  /."A-CC*  br.Mf  V 

15AO  PHIM  •  Fl*.^'  Fl/i:cC'\»' Fl/SLCtb* 

1550  Fph  K=  1  10  33 

15fO  PklN'l  v.CK>,V(K  ^hC}<)^CCK.‘ 

157J  kiiAl  K 

1700  LAlA  1 .  1  7b>  .  49  7»  .  <’b>9o 
^710  OaIa  -bb^.b^D 
*  "''bo  LhIa  loo 
9999  c'\‘o 
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